Continuous Spectra of Generalized Kronig-Penney Model 
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The standard Kronig-Penney model with periodic 5 potentials is extended to the cases with 
generalized contact interactions. The eigen equation which determines the dispersion relation for 
one-dimensional periodic array of the generalized contact interactions is deduced with the transfer 
matrix formalism. Numerical results are presented which reveal unexpected band spectra with 
broader band gap in higher energy region for generic model with generalized contact interaction. 
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I. INTRODUCTION 

The contact interaction occupies a special position in 
quantum mechanics. The system with contact interac- 
tion is often rigorously solvable and also is useful in 
examining the effect of small obstacle on particle mo- 
tion. The first influential work on contact interactions 
was done by Kronig and Penney The model, which 
has potential consisting of a periodic array of S functions, 
has been widely regarded as a standard reference model 
in the solid-state physics for more than half a century. 

In spite of the seeming simplicity of contact inter- 
actions, there are several non-trivial aspects which are 
largely left unexplored. Even in the simplest setting of 
one dimension, there has been a historically longstanding 
problem of realizing the generalized contact interaction 
in the small-size limit of a local self-adjoint interaction. 
In one dimension, there are a four-parameter family of 
generalized contact interactions which conserve the cur- 
rent at both sides of the interaction |3( • This corresponds 
to the fact that the one-dimensional kinetic energy oper- 
ator T = --^i with domain Cq°(R\{0}) has deficiency- 
indices (2, 2) 2 Q. In non-relativistic formalism, the cur- 
rent operator is given by 



j(x) = &(x)a 2 i&(x), 



(1) 



where ^ is defined in terms of the wave function tp and 
its space derivative as 



(2) 



The matrix 02 in Eq. Q is the second component of Pauli 
matrices; 
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In Eq.©, m is the mass of a particle. If we put a contact 
interaction at x = 0, the connection of between both 
sides of the origin can be characterized by 



where the current conservation j(+0) = j{— 0) requires 

VW = <7a. (5) 
The generic solution of Eq. J5J is given by 
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where 6 € R and — (38 = 1. The condition © covers 
a four-parameter family. (See 0, 0, E| f° r the discussions 
on rigorous characterization of the full parameter space.) 
In addition to the usual 5 potential 
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this family includes the so-called e potential (also known 
by a misnomer, 5 potential). The e potential induces 
such a boundary condition that the wave function has 
continuous first derivative on the right and left, but it 
has a jump proportional to the first derivative [jj; 
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(8) 



*(+0) =V*(-0), 



(4) 



In we have constructed the e potential in the small 
distance limit of three nearby 6 potentials and realized 
in terms of usual 5 and e potentials a three-parameter 
family of self-adjoint extensions under the assumption of 
time reversal symmetry. This treatment have been gener- 
alized with finite range potential with double short-range 
limit with different scales [jj. Also, the relativistic origin 
of e potential has been discussed in ^FJ- Contrary to 
the general belief, some of the non-<5 contact interactions 
appears in realizable settings in quite a natural manner 
. and, in fact, they are essential in understanding the 
nature of one-dimensional many-body system |l2j . 

In light of these facts, it is quite interesting to examine 
the nature of one-dimensional system with a periodic ar- 
ray of generalized contact interactions. For this purpose, 
we generalize the Kronig-Penney model to the cases of 
generalized point interactions. We are particularly inter- 
ested in the question whether the original Kronig-Penney 
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FIG. 1: Periodic potential of Kronig-Penny model with gen- 
eralized contact interaction. The two-by-two matrix V £ 
SL(2, R) represents the connection condition at each inter- 
action in the transfer matrix formalism. 



model with ^-interaction is a generic zero-range limit of 
more realistic models with finite-range interactions. In 
this paper, we assume that the system has time-reversal 
symmetry. In this case, we have 9 = in Eq.©, hence 
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After giving the formulation of the generalized Kronig- 
Penney model with generic contact interactions in Sect. 2, 
some elementary numerical results are presented in 
Sect. 3. It will be shown that the standard periodic array 
of 5 potentials has a specific band structure, compared 
to the generic cases; The band gap tends to disappear in 
the high energy limit for periodic 5 array, while the band 
width becomes narrower for higher bands in generic cases. 
The current work is summarized in Sect. 4. 



II. FORMALISM 

We consider a one-dimensional periodic array of a gen- 
eralized contact interaction, the connection condition of 
which is described by V in Eq.©. We assume that the 
interactions are located at x — na, (n = 0, ±1, ±2, • • • ). 
Here we denote the lattice interval by a. The assumed 
potential is shown in Fig.l. 
Schrodinger equation is given by 
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<p"{x) = Eif(x), (x =/= na) 



(10) 



with the boundary condition 

*(na + 0) = V*(na-0) (11) 

at x = na. In the transfer matrix formalism, Ea. HlOjl is 
written as the first-order coupled equation; 



*'0) = H9(x), {x^na) 
where we define 
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(12) 



(13) 



with ko = \/2mE. The solution of Eq. I|12|) is written as 
*(a:) = G{x-x Q )*(x ) (14) 
by using the exponential function of Tlx; 



Q(x) 
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cos(fc a;)/2 + 



sin(fcoT) 
k 



n. (15) 



Here we assume na ^ [xq, x] in Ea. l(T4*|) . In Ea. (|T3l) . I2 is 

the two- by-two identity matrix. Note that det Q (x) = 1 
because of TrTC = 0. 
We can see 



G(x)u ±kQ = e ±lk ° x u ±ko 



with 



u ±fc 



±( 1 



(16) 



(17) 



Namely, Q{x) has eigenvalues e ±lk ° x with the associated 
eigenfunction u±fc in Ea. (|17fl . Also the complex conju- 
gate of 0(x) satisfies 



with 



G\x)v ±ko =e* ikax v ±ko 



v ±fc0 - ^ ( ± 



(18) 



(19) 



That is, Q^(x) has eigenvalues e Tlk ° x with the associated 
eigenfunction v± k(l in Ea. H19(l . The eigenfunctions u± ko 
and v± ko satisfy the bi-orthogonal relations; 

vL u ±fc = 1, v^ fco u ±fco = 0. (20) 
By Bloch theorem, we can set 

cp(x) = e ikx u(x), (k G R) (21) 
in Ea. HlOfl . where the function u has period a; 

u(x + a) = u(x). (22) 

Since 

v '(x) = e lkx (u'(x)+iku(x)), (23) 
cp"( x ) = e ikx (u"(x) + 2iku'(x) - k 2 u(x)), (24) 

we obtain the equation for u; 

- u"(x) - 2iku'(x) + k 2 u(x) = k%u(x), (25) 

(x 7^ na). In the vector notation 



*(*) = ( v 



we rewrite Eq. (J25J) as 

&(x)=H9(x), (x^na) 



(26) 



(27) 
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with 



n 



2m 



-2ik 



(28) 



It is clear from Eas. (|21|) and (I23|) that the two vectors 
*&(x) and *&(x) are related by 



*(x) = M^>{x) 



with 



M = e lkx 



1 
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The solution of Eq. ***7|) is given by 

#(a>) =g(x-x )$(x ), 
where Q(x) is the exponential function of Tlx; 



g(x) 



(29) 



(30) 



(31) 



(32) 
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Here we assume na (jt [xq, x]. Since TrTL = —2ik, we have 



detQ(x) 



D (TrH)x 



-2ikx 



It is seen from Eas. pi(l 



and Ij29(l that the connection condition for \P at x = na 
is given by 



with 



*(na + 0) = V*(na-0) 



ifc i J ^ ( _ik_ i 
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(33) 
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Note that detV = detV = 1. The periodicity for u in 
Ea. (|2*2l is equivalent to 



*(x + a) = 
In particular, with x — na — 0, we have 

£?(a)V*(wz - 0) = *(na - 0), 



since 



#((n + l)a-0) = g~(a)ty(na + Q) 

= g{a)VV{na~ 0). 



(35) 



(36) 



(37) 



Ea. (|36fl requires the two-by-two matrix Q{a)V has eigen- 
value 1: 



det(/ 2 - G(a)V) = 0. 



(38) 



Since det(<?(a)V) = det^(a) = e ~ 2tka , the other eigen- 
value of the matrix g(a) V should be e~ 2lka . This indi- 
cates 

Tr{g(a)V) = l + e~ 2lka . (39) 
A simple matrix calculation shows 

MUAr 1 = -ikh + U, (40) 
which gives the relation between g{x) and g(x); 



Mg{x)M~ 1 = 



^MHM-^x _ e -ikx 



c g(x). (4i) 



From Ea. (|4TT> . we obtain 

g(a)V = e'^M^g^VM. 



(42) 



Thus we conclude that the condition P9")l is equivalent 
to 



Tr(g(a)V) = 2cos(ka). 



(43) 



Eq. (|43|l determines the dispersion relation for a periodic 
array of the generalized contact interaction characterized 
by the connection condition 0. Inserting Eqs.@ and 
(|TK|l into Eg. 1(4*31) . we obtain the eigenvalue equation 



(a + 7) cos(fcoa) + sin(fcoa) ( -, — P — 77-^ 

1 Kq 2m 



2cos(ka). (44) 



III. NUMERICAL EXAMPLES 

In this section, we give numerical examples of the band 
spectrum in several cases; the usual periodic <5, periodic e 
as well as periodic array with some typical one-parameter 
families of generic contact interactions. It will be shown 
that the band structure of periodic S potential is not 
generic; The band width becomes broader even in the k 
(wave number) space for 6 array as the energy increases, 
contrary to the generic cases. Throughout this section, 
we take the mass of particle m = 1/2 and the lattice 
interval a = 1 in numerical calculations. 

Let us begin with the periodic 6 potentials. From 
Eq . (|44() together with Eq.(|7j), we obtain the eigenvalue 
equation for S array; 



cos(fco) 



v sin(fco) 

2kn 



cos(fc), 



(45) 



where v is the strength of each S. Ea. i|45|) is the 
well-known result in the standard Kronig-Penney model. 
Fig. 2 (a) shows the band spectrum as a function of v. 
With the parameters v and E — k^ in the shaded region, 
the equation l|45|) has the solution k. The condition for 
the existence of the solution is given by 



cos(fc ) 



v sin(fco) 



2fcn 



< 1. 



(46) 
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FIG. 2: Band spectrum for a periodic array of generalized 
contact interaction characterized by one-parameter family of 
connection conditions; Eqs.(7), (8), (49), and (51) in (a), (b), 
(c), and (d), respectively. In all cases, the horizontal axis is 
the parameter associated with each family, while the vertical 
fe 2 

axis is the energy E = w 2 -. The mass is set to m = 1/2. The 
lattice interval is set to a = 1. 



For v — 0, one sees the continuous spectrum for positive 
energy and there is no bound states, as expected. 

For v 0, the band structure appears and the band 
width becomes narrower as the strength \v\ increases. In 
the limit of \v\ — ► oo, one obtains point spectrum E — 
(nir) 2 . For negative v, there is a negative energy band 
which goes to the minus infinity in the limit of v — ► — oo. 
In both limits v — > ±oo, the boundary condition around 
each S corresponds to the so-called separated boundary 
condition. Each region between two neighboring <5's are 
separated from the other regions and the wave function 
satisfies the Dirichlet boundary condition such that it 
vanishes at each <5. 

For e potential of strength u, we obtain the eigenvalue 
equation 



cos(fco) 



uko sin(fco) 



cos(fc). 



(47) 



Fig. 2(b) shows the band spectrum as a function of u. The 
condition for the existence of the solution is given by 



cos(fco 



uk sin(fc ) 



< 1. 



(48) 



The case of u — corresponds to the free space. For 
u 7^ 0, the band structure appears and the band width 
becomes rapidly narrower as the strength \u\ increases. 
In the limit of |u| — > oo, one obtains point spectrum E = 
(nir) 2 . In this limit, each region between two neighboring 
e's are separated from the other regions and the wave 
function satisfies the Neumann boundary condition such 
that its derivative vanishes at each e. For negative u, 
there is a negative energy band which goes to the minus 
infinity in the limit of u — ► —0. 

To see the generic cases, we show in Fig2.(c) the band 
spectrum for the one-parameter family of the connection 
condition: 



V = 



cosp — sinp 
sinp cosp 



( — 7T < p < 71"). 



(49) 



The band spectrum is determined by 



cos(p) cos(fc ) 



sin(p) sin(fco) f 7 . 1 

— 2 — r + fcu 



< 1. (50) 



With p 0, 7r, we see the band structure as in other 
generic cases. The band spectrum has period tt, as seen 
from Eg. (|50|) and the energy spectrum for p = tt is the 
same as in the free space (p = 0). However, the wave 
function for p = n differs from the free one. Indeed, it 
has the same amplitude as the continuous wave function 
in the free space, but it changes the phase and as a result 
it is discontinuous at each obstacle. Such "duality" of 
p = and p = 7T induces double spiral structure in energy 
spectrum [l3j, as implied in Fig. 2(c). 

Fig. 2(d) shows the case for the one-parameter family 

of 



y f coshp sinhp \ 
I sinhp coshp J 



( P e R). 



(51) 
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In this case, the band spectrum is determined by 
sinh(p) sin(fc 



cosh(p) cos(fco) 



k - 



< 1. (52) 



As in other cases, the band width becomes gradually nar- 
rower, as the perturbation becomes larger. 

An important point is that the S potential, among 
generic point interactions, has rather a special energy 
dependence of band structure. To see this, we show 
the band spectrum including higher energy region in 
Figs.3(a)-(d), which correspond to Figs.2(a)-(d), respec- 
tively. In Fig. 3, the vertical axis is the wave number ko, 
instead of the energy E, which is suitable for our present 
purpose, since the state density is constant [of 0(k )} in 
one dimension. One can see that in generic cases includ- 
ing e, the band width becomes narrower for higher bands, 
while the band gap tends to disappear for S array. 

The reason is easy to understand by examining the scat- 
tering properties by a single contact interaction. 

In the transfer matrix formulation, the bi-orthogonal 
eigenvectors (|17fl and l|19|) serve to examine the scatter- 
ing properties. Since e ±ik ° x u±k are the solutions of the 
equation H12|) . the wave function is written as 



T s e ik ° x u +ko , 



R S e 



— ikox 



u -feo 



(x < 0), 
(x > 0), 



(53) 



where Tg and Rs are transmission and reflection coeffi- 
cients, respectively. We assume that the contact interac- 
tion is placed at x = and the incident wave comes from 
minus infinity in Ea. l|53[) . From the connection condition 
(@J, we obtain 



Tu +ko = V(u +feo + i? S U_ 



(54) 



Multiplying v^ fc V 1 from the left and using the bi- 
orthogonal relations (|20|l , we can estimate the transition 
probability as 



\T\ 



U +fer 



a 2 + Y + 2 + 6 



k 2 
4m 2 



P 2 



Am 2 



k 2 



The reflection probability is given by 
|i?| 2 = 1- \T\ 2 . 



(55) 



(56) 



Ea. (|55|l shows \T\ 2 — > as ko — > +oo, namely perfect 
reflection in the high energy limit for generic cases. The 
exception arises in case of S = 0, namely S potential. In 
this case, we have 



|r a | a 



2 21 - 1 

v m 



k 2 
k 2 



(57) 
(58) 
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FIG. 3: The energy dependence of the band spectrum. The 
indications are the same as in Fig. 2 except that the vertical 
axis is the wave number ko- 
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Thus, \Tg\ 2 — > 1 as fco — > +00, namely perfect trans- 
mission is realized. This explains why the band width 
broadens even in the fco (wave number) space for the 5 
array as the energy increases and the band gap disap- 
pears in the high energy limit. 

In the low energy limit, the e potential rather than 
5 shows a special nature. From Eas. (|55|l and (|56|l . one 
recognizes in generic cases, |T| 2 — > as fco — ► 0, i.e. 
perfect reflection in the low energy limit. The exception 
arises in case of /3 — 0, namely e potential. In this case, 
we have 



\T K 



u 2 k 2 



16m 2 
16m 2 



1t 2 fcg 



(59) 
(60) 



Thus, |T £ | 2 _ — 
sion is realized. 



1 as fco 



0, namely perfect transmis- 



interactions in one dimension. The transfer matrix for- 
malism serves to deduce the eigenvalue equation which 
determines the dispersion relation. 



Numerical analysis of the band spectra shows that the 
band structure by the usual periodic 5 array is not generic 
and the band width becomes broader in the fc space as the 
energy increases, while it tends to be narrower in generic 
cases. This fact opens up an interesting possibility that 
several distinct classes of band structures found in var- 
ious materials may be modeled and classified in terms 
of simple and solvable but versatile model of general- 
ized Kronig-Penney hamiltonian. It would be of interest 
to extend the current approach to related problems with 
direct experimental relevance, such as the Wannier-Stark 
ladder problem 0, 0] . 



IV. CONCLUSION 

We have formulated the generalized Kronig-Penney 
model which has a periodic array of generalized contact 
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